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^1/ ' Abstract: Chiral condensate and r]' meson mass spectrum are studied under the influence 

^SJ ■ of an external Abelian magnetic field. We work within the D3/D7 Karch — Katz model of 

f~^ ■ flavoured AdS/CFT with supersymmetry broken by the Constable — Myers deformation 

^ . of the metric. It is shown that this setting yields an analytic (quadratic) dependence of 

f^ , condensate on fleld, typical for the Nambu — Jona-Lasinio model, rather than the non- 

analytic (linear in field) result, typical for chiral perturbation theory in the exact chiral 
limit. We conjecture that the analytic (quadratic) result might be put into correspondence 



\-{ ' with the leading-order in the 1/Nc decomposition for the condensate. This leading order 

in the 1/iVc approximation has not yet been derived from the chiral perturbation theory. 
Thus the dual model yields the quadratic field dependence of the condensate, which is 
beyond the range of feasibility of chiral perturbation theory. 
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Introduction 

The behaviour of QCD vacuum in strong electromagnetic fields has recently attracted a 
great deal of attention (e.g. |||, |2|), reinvigorating the subject which had been started byQ. 
Lattice simulations |Q, P, Simonov's string model Q are just a few of the recent studies 
of QCD vacuum in external fields to be mentioned here. In this work we try to describe 
the behaviour of mesonic spectra and condensate from the perspective of duality. This 
article is organized as follows. In the following Section |^ treatment of mesonic masses' and 
condensates' in external fields is reviewed. It is explained why traditional field-theoretical 
approaches, are still demanding a non-perturbative insight, possibly coming from the realm 
of dual models. Then in Section y a short description of the specific dual model is given, 
which we are going to apply. In the subsequent Section ^ the numerical calculations are 
presented. We conclude in ^. 

1. Motivation 

The QCD vacuum is quantitatively described by its chiral condensate, gluonic conden- 
sate, pion decay constant and some other physical quantities. Below we shall revisit the 
properties of some of these objects in the electromagnetic background. 



1.1 Condensate 

QCD chiral condensate {qq) is the order parameter of chiral symmetry breaking. Important 
ideas of chiral symmetry breaking catalysis were being developed by Gusynin, Miransky 
and Shovkovy. In |2|, ^, |8| they study enhancement of chiral condensates in 2 + 1 and 
3 + 1-dimensional Nambu — Jona-Lasinio-like models. 

The issue of condensates in an external magnetic field was resolved by Schramm, 
Mueller and Schramm H, and by Smilga and Shushpanov in [10|. For small magnetic 
fields H, and in an exact chiral limit 

, , f eHln2\ , , 

(wk = (w)o(^i + Y^^J- (1-1) 

Note that the linear term in H has a -^ factor, for /^ ~ A'c A second-loop correction to 



this result was calculated by Shushpanov and Agasian |]ll|. It is instructive to compare 



this low-energy QCD computation with a Nambu — Jona-Lasinio model computation made 



by Klevansky and Lemmer [12| 



(9-g)H = (g-5)o(l + c^^|^), (1.2) 



where c is some model-dependent coefficient. The linear dependence ( |l.lD by Smilga and 
Shushpanov is non-analytic (has a square-root type cut) in terms of the invariants of the 
external field, i.e. is organized as ~ yF^ . This might seem to be inconsistent from a 
first view. However, this non-analyticity is of vital importance. It means there are no 
other massive parameters in the low-energy domain, where the chiral perturbation theory 



is valid. The non-analyticity of (1.1) is a direct signature of vr- meson being a Goldstone 
particle. If chiral limit is violated, the dependence will be analytic. One must work here 
in the exact chiral limit, for otherwise all other massive hadronic states in the vacuum 
energy loops must be taken into account. Nambu — Jona-Lasinio model is at the same time 
seen to be deficient to describe full QCD, as it does not reproduce the correct non-analytic 
behaviour of the condensate, representing the Goldstone particles. 

Chiral condensates in arbitrary electromagnetic fields were calculated by Cohen, Mc- 



Gady and Werbos in |13|. They have obtained expressions for electric, magnetic, and 
arbitrary configuration of constant fields. Their results are basically obtained in the same 
Heisenberg-Euler technique type as those of Smilga and Shushpanov, and perfectly repro- 
duce the latter as a particular case. 

1.2 Limitations of traditional approaches 

The above chiral perturbation theory results have the status of exact low-energy theorems. 
However, they have their domain of applicability, as explained in the review paper by 
loffe ||l3|. The one-loop result has been reminded above. This means there will be next- 
order loop corrections in chiral perturbation theory to this value. Chiral perturbation 
theory has also limitations due to the fact that quarks' and gluons' degrees of freedom 
are fully absent in it. Therefore, other models have to be considered to be compared with 
chiral perturbation theory estimates. 



A class of modern (supersymmetric and non-supersymmetric) QCD models, which 
would be natural to test for the behaviour of condensates and meson masses, are AdS/CFT 
models with flavours. For a review on AdS/CFT with flavours or in non-supersymmetric 
backgrounds see [l^, 16|. They are generally constructed basing on Maldacena's conjec- 



ture 1 17] on equivalence of the IIB type closed string theory in the bulk of AdS^ xS^ and 
AA = 4 SYM theory in the four-dimensional flat spacetime. For a review of Maldacena 
conjecture in general, see [^, ^. QCD is, of course, a non-supersymmetric and a non- 
conformal theory, so various symmetry breaking techniques are applied to make the model 
resemble the reality. The two most common approaches are: "bottom-up" approach, and 



"top-down" approach. A "bottom-up" construction, see e.g. |2^ is usually constructed 
with a 5-dimensional action, which one has to "guess" , so that it fits as many QCD results 
as possible. On the other hand, the "top-down" approach (e.g. ||2l|]) is constructed start- 
ing with a 10-dimensional geometrical setting, in which special elements are supposed to 
reproduce the dynamics of QCD degrees of freedom. 

Flavoured AdS/CFT correspondence in an external magnetic field was studied by 
Filev et al. in |^2[ . They produce a spectrum of mesons from pure-AdS background, which 
satisfies the Gell-Mann — Oakes — Renner relation. In |2^] thermodynamic properties of the 
gauge theory in a magnetic field have been studied in the same framework. Properties of 



the theory in an electric field were obtained in |24] by the same method. 



AdS/CFT with flavours in external fields and at finite temperatures have also been 
studied in p5|. The authors calculate a number of external-field-dependent properties for a 
supersymmetric background, such as meson masses in electric and magnetic fields. Sakai — 
Sugimoto model in external fields was studied in ^^. It has been concluded that Sakai — 
Sugimoto model is consistent with the picture of magnetic catalysis of chiral symmetry 
breaking. Phase transitions in Sakai — Sugimoto models due to switching on of electric 
and magnetic fields were discussed in [ES]. Pair production in an electric field in Sakai — 



Sugimoto model was studied in |3C]. When the present paper was being completed, two 
works on holographic QCD at finite temperature, magnetic field and chemical potential 
appeared on the same day ||3^, 32]. This is the evidence for the great interest to the 



different aspects of AdS/CFT in external Abelian fields that is present nowadays. 

2. D3/D7 model in Constable — Myers background with a Kalb — Ramond 
field 

In this short letter a very simple model is discussed, which features many of the basic QCD 
characteristics: confinement, conformal symmetry breaking and spontaneous chiral sym- 
metry breaking. We put it into a magnetic field and observe the behaviour of condensates 
and mass spectra. 

Below we follow what is known as Karch-Katz model with Constable — Myers defor- 
mation. We rely in this passage essentially (sometimes literally) on [^. This geometry 
conjecturally describes a AA = 4 SYM broken by non-zero expectation values for all S0{6) 
singlet operators. It also inherits from pure Karch — Katz model a pack of D7 probe branes, 
which do not affect the metrics. This requires us to work in the quenched approximation 



Nf ^ Nc- The Constable — Myers metric is organized as 
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^= il + TTl -1- (2.2) 



This form of the metric makes it easy to see that it behaves asymptotically at r ^ oo as 
pure AdS, but differs from it near the singularity. 

The Constable — Myers solution requires a non-trivial dilaton as well 

^20 ^ ^200 ("l + ^"j (2.3) 



r^ 



and a C4 form field 



C(^4^) = --H ^dxo Adxi Adx2 /\dx3, (2.4) 



4 
with conditions imposed upon deformation parameters 

A'^ + 6'^ = 10, 
^=2F- 



(2.5) 



For a more convenient embedding of the D7 brane, a coordinate transformation is per- 
formed, which will explicitly separate the 4-dimensional and 6-dimensional hyperplanes: 



(2.6) 
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and the dilaton is 








U4_54) 







(2.7) 



(2.8) 

We have already learned about the equivalence of Kalb — Ramond field in the bulk and the 
Maxwell field on the brane. The D7 brane does not change the metric in the quenched 
approximation. The dynamics of the brane is described by a Dirac — Born — Infeld action 



d^^Jdet (27rB^p + Ina'F^^ + g^,—-—^\ + I d^^C^ AFAB (2.9) 

Here B^^ is the Kalb — Ramond field, defined in the bulk, which is projected to the brane as 
Bai3 and Fap is the usual Maxwell field on the brane. A constant field -F23 = —-^32 = -B is 
chosen, all other field components being zero. The Chern — Simons term does not influence 



classical dynamics. It may give a contribution into the oscillations describing mesonic 
masses. Further the embedding geometry and gauge condition are specified. The D7 brane 
runs through the directions of coordinates xo,xi,X2,X3,wi,W2,ws,W4. These coordinates 
are respectively ^i . . . ^^s internal coordinates of the brane world- volume. It doesn't run 
through the remaining W5,wq. The latter coordinates are embedding coordinates of the 
brane into the targetspace. They are functions of .^j. Solutions in the form w^ = w{p),WQ = 
will be sought, where 



'wj + wl + wl + wj. (2.10) 

With such an Ansatz and in the metrics given above, the DBI action is organized as 

S = -fi7Jd^^G{p,w)^/l + w'^ip)^Jl + Bygf,, (2.11) 

where 

The equations of motion will look like 

"^ ^ ^'^' 'l + SVfffi) - ^^+^^ (g^I + Bygl^ = (2.13) 



dp \y/l + w'"^ 

They are solved them numerically in the next section. First quark masses and condensates 
are extracted and fitted with appropriate interpolation functions. Then small oscillations 
around the classical solutions are studied. The spectra of these oscillations are identified 



with meson masses according to known rules [27| of AdS/CFT correspondence. 



3. Condensate and spectra 

3.1 Condensate 

The standard lore is: one must search for physical solutions of these non- linear second- 
order differential equations, which have the following asymptotics in the infinity: 

w(p) =m+ ^T. (3.1) 

P 

Then the parameters m and c correspond to quark mass and chiral condensate: 



m — "^ 



{m) 



c 

(27ro')^ ' 



(3.2) 



where a' is string tension parameter. Contrary to the physical solutions, unphysical ones 
are those ending in the singularity of Constable — Myers metrics, or going to infinity at 
p — > 0. The singularity is marked by an ellipse denoted "singularity" in Fig. (||). Physical 
solutions can be defined by boundary condition w'[{)) = 0. It happens that the generic 
solutions are unphysical ones. 
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Figure 1: Dependence of condensate on magnetic field and mass. 
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Figure 2: Different embeddings of the spectator D7 branc into Constable — Myers background. 



To obtain physical solutions, one imposes 



w'{0) = 0, 

w{0) = 11)0 = const. 



(3.3) 
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Figure 3: Magnetic catalysis of chiral symmetry breaking in Karch — Katz model with Constable- 
Myers deformation, exact chiral limit m = 0. 
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Figure 4: Mass spectra of rn^ as functions of B. Thin lines on the left (m = 1) and right 
(771 — 0.1) plots show interpolation Jm^ = aB and Sm'^ = aB^. One can see that neither of these 
interpolations is satisfactory. 



For each value of wq above some value Wmin ~ 1-34 this will yield a curve Fig. (|2|), asymp- 
totic behaviour of which will reveal some definite m and c. This allows one to build the 



dependence of condensate on quark mass Fig. (|l|). Doing the same thing with different val- 
ues of magnetic field B one gets a shifted curve. It's a subtlety of this method that in order 
to understand how the condensate shifts in the field, one must take a section of Fig. (|l]) at 
a constant m rather than follow some definite point in the plot. The resulting dependence 
of condensate on the field is shown in Fig. (|3|). It is noted here that the results for B = 
coincide with those calculated in the same background in |3S]. 



We analyze the "experimental" dependence. It is natural to expect either linear (as in 
true QCD) or quadratic (as in NJL) condensate growth with the field. In our case, approx- 
imation with a quadratic polynomial comes out to be quite effective. In the picture Fig. (g|) 
one can see the comparison between the linear and quadratic approximations, and judge 
in favour of the latter. 

This quadratic dependence on field value corresponds very nicely to the picture of 
magnetic catalysis of chiral symmetry breaking in [^, ||] in NJL models. On the other 
hand, it does not correspond to linear condensate shift, predicted by the low-energy QCD 
effective action by Smilga and Shushpanov [10|. This phenomenon may be given a nice 



qualitative explanation. Condensate expression ( |1.1| ) is a part of the series in powers of 
for /^ ~ y/N^. It starts with the jt- term. There may be a term, dependent on field, 



and containing -^ in the zeroth power. To our best knowledge, such terms have not been 
reported in chiral perturbation theory. On the contrary, dual models restore the missing 
leading-order -^ contribution. 

3.2 Meson spectra 

Small fluctuations of the classical solutions to the equations of motion govern mesonic 
spectra. There are two types of these fluctuations: those corresponding to Goldstone 
(in the large Nc limit) mesons eta' , and those corresponding to the non-Goldstone ones. 
The former are fluctuations of the angular coordinate in OxgXg plane, the latter are the 
fluctuations of the radial coordinate. The equations for Goldstone part of the spectra 
are [p3l: 



A 

dp 



G^l+Bygl, ' ^^ ■ ' " ■■ ' ' 



Vl+w''^ P 



dpfip) 






This is a Sturm — Liouville eigenvalue problem on function f{p), which must be solved with 
the following boundary conditions: 

/•^'<°' = » (3.5) 



1 



/(p)lp^oo -^ -71 



P 



The function w in the equation ( p.4| ) must be taken from the previous section. These 
equations are solved in Mathematica environment by shooting method. One starts with 
solutions behaving like -^ in the infinity, and step-by-step find the value of M, which yields 
the desired behaviour at the origin. Again, the resulting regular solution is never reached, 
but it can be approximated as a separatrix of solutions singular at origins, to any desired 



accuracy. In our calculations we have obtained M^ with 4 decimal digits. The results for 
Goldstones are shown in Fig. (Q), where 

5ml = ml{B)-ml{0). (3.6) 

We have tried, basing on previous experience both on field theory side and gravity side, to 
approximate the field dependence of 6m'^ by either linear drn^ ~ i? or quadratic Sm'^ ~ 
B^ dependence. However, our numerical analysis has shown that neither can be a valid 
approximation, even for small B, which is shown in Fig. @. Comparing this to the linear 
dependence for masses in the chiral limit of chiral perturbation theory and to the quadratic 
dependence given for pure AdS background in ||2^, we conclude that dynamics of our model 
might be away from both the predictions of chiral perturbation theory and pure AdS model. 

4. Conclusion 

A qualitative conclusion can be drawn upon analyzing the dependence of condensates on 
the field. We can see that the linear field dependence of QCD condensate from chiral 
perturbation theory is not reproduced at all. Instead, a quadratic dependence is retrieved. 
Our conjecture to explain this phenomenon is very simple. Chiral perturbation theory 
estimate, as given in the cited references, misses the leading-order in ;^. It starts with 
the next-to-leading order in -^. On the other hand, duality might reproduce the leading- 
order effect. Nevertheless, the search for a true dual model of QCD must still be in 
progress, for meson mass spectra cannot be easily given a qualitative explanation. One of 
possible improvements of the model would be to take into account back-reaction effects. 
In our setting, the D7 brane was a probe brane, a self-consistent supergravity solution in 
a background of a stack of D3 branes and a D7 brane would be more complicated. 
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